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Vibrationally coupled electron transport through single- 
molecule junctions is considered. Reviewing our recent 
theoretical work, we show that electron-hole pair cre- 
ation processes represent the key to understand the vi- 
brational excitation characteristic of a single-molecule 
contact. Moreover, these processes can lead to a number 
of interesting transport phenomena such as, for example, 
negative differential resistance, rectification, mode-se- 



1 Introduction The study of single-molecule junc- 
tions represents a vibrant and challenging field of research 
comprising fundamentally important aspects of quantum 
and many-body physics as well as the possibility for tech- 
nological applications in nanoelectronic devices f 1 1. Elec- 
tron transport through a single-molecule junction can be 
understood in a similar way as through quantum dots, in- 
cluding electron-electron interactions [2 1 and/or quantum 
interference effects 014]. In addition, the vibrational de- 
grees of freedom of a molecule play an important role |]5] 
[6]|4], including, in particular, nonequilibrium effects such 
as current-induced vibrational excitation [7 8 9. T0][TnfT2l 
[T3l[T4 15 164 17|. This is in contrast to quantum dot sys- 
tems fl8l . where phonon degrees of freedom are often 
considered as being part of the environment rather than as 



lective vibrational excitation and a pronounced tempera- 
ture dependence of the electrical current. Thus, electron- 
hole pair creation processes are crucial to elucidate 
the basic mechanisms of vibrationally coupled electron 
transport through a single-molecule contact, despite the 
fact that these processes do not directly contribute to the 
electrical current that is flowing through the junction. 
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active degrees of freedom, and is associated to the small 
size and mass of a molecular conductor. 

The key to understand this nonequilibrium transport 
problem is to account for all relevant processes. Quite intu- 
itively, this includes transport processes, where an electron 
is transferred from one electrode to another. In the pres- 
ence of an external bias voltage, these processes give rise to 
an electrical current that is flowing through the molecule. 
However, these are not the only processes that occur in 
single-molecule junctions. Considering the simpler setup 
of a molecule adsorbed on a surface, it is well known that 
electron-hole pair creation processes, where an electron 
transfers from the substrate onto the molecule and back, 
constitute the most important relaxation mechanism for vi- 
brational excitation |[T9| . Such pair creation processes do 
also occur in single-molecule junctions, although their im- 



Copyright line will be provided by the publisher 



R. Hartle, U. Peskin and M. Thoss: Electron transport through molecular junctions: Role of pair creation processes 



portance for the respective nonequilibrium transport prop- 
erties has been established only recently ||9) |20]|2T] ll6,22, 
[T7P1l. In this article, we review the most important as- 
pects of our work on this subject i20]|2T1[T6]|22l|23][T7]|4l . 
This includes a detailed understanding of the vibrational 
excitation characteristic, in particular the counter-intuitive 
phenomenon that larger levels of vibrational excitation are 
obtained for systems with weaker electronic -vibrational 
coupling 1141221 . In addition, pair creation processes may 
also have a substantial influence on the electrical trans- 
port properties of a single-molecule contact, giving rise, 
for example, to negative differential resistance, rectifica- 
tion |[T4l and important implications for spectroscopical 
applications of single-molecule junctions [23 1. Pair cre- 
ation processes can even facilitate a mechanism to con- 
trol the excitation level of specific vibrational modes se- 
lectively via the external bias voltage |20)2T] and, in the 
presence of destructive interference effects, of the electri- 
cal current via the temperature in the electrodes H22JI171 . 
Note that the latter has been experimentally verified only 
recently J4). 

The article is organized as follows. In Sec. [2] we in- 
troduce the model Hamiltonian and the nonequilibrium 
Green's function (NEGF) approach ll24ll25lll0ll2Tl that 
we use to describe electron transport through a single- 
molecule junction. Our results are presented in Sec. [3] 
where we set the stage by a discussion of the basic trans- 
port processes (Sec. 13. 1) and electron-hole pair creation 
processes (Sec. 13.2b . including a detailed analysis of the 
vibrational excitation characteristic of a generic molecu- 
lar junction. Transport phenomena that are attributed to 
the influence of electron-hole pair creation processes are 
summarized in the last section, Sec. 13.31 Note that the 
presented data has already been published in Ref . 11231 . 

2 Theory 

2.1 Model Hamiltonian We consider a model for a 
molecular junction obtained by a partitioning of the over- 
all system into the molecule and the left and right leads. 
The molecular part of a single-molecule junction (M) is 
often described by a set of discrete electronic states ||7][T] . 
The Hamiltonian of the corresponding electronic degrees 
of freedom can be written as Ifl4l 



Hei = X^ e m cl 



(1) 



mEM 



/ UmnxCrnCm OjyijyC^Cji u n J, 



m<n£M 

where e m denotes the energy of the mth molecular state 
and c^lcm the corresponding creation/annihilation oper- 
ators. Thereby, the index m distinguishes, in principle, 
different molecular orbitals, including the spin of the elec- 
trons. However, as the effects and mechanisms that are 
discussed in this review article do not explicitly involve 
spin degrees of freedom, we suppress them in the fol- 
lowing. Coulomb interactions are described in the model 



Hamiltonian H by Hubbard-like electron-electron interac- 
tion terms, C/„ m (cJ„c m - 5 m )(clc n - S n ). Thereby, we 
account for the fact that the single-particle energies e m 
are determined with respect to a specific reference state 
of the molecule (e.g., the ground state of the uncharged 
molecule) II26II27I and that, effectively, these energies 
include Coulomb interactions between the electrons of 
the reference system. Therefore, we distinguish occupied 
(6 m = 1) and unoccupied states (5 m = 0) of the refer- 
ence state to avoid double counting of electron-electron 
interactions. 

Besides the electronic degrees of freedom, nuclear mo- 
tion and/or electronic-vibrational coupling plays a key role 
in transport through molecular conductors [7, 10l28l l5]l4l. 
which is associated with their small size and mass. We take 
into account the vibrational degrees of freedom of a molec- 
ular conductor as harmonic oscillators that are linearly cou- 



pled to the electron (or hole) densities (c 



$m) 



271 



#vib — 2^ ^a«i a « 



) ^mflVftl'-m^ Om)j \^-J 



where the operator aj, denotes the creation operator of the 
ath oscillator with frequency Q a and Q a = a a + a' a the 
corresponding vibrational displacement operators. The re- 
spective coupling strengths are denoted by X ma . Note that 
we identify the vibrational degrees of freedom of the junc- 
tion as the normal modes of the aforementioned reference 
state. As a consequence, electronic-vibrational coupling is 
required to vanish for that state, which, similarly as for the 
electron-electron interaction terms, is ensured by the pa- 
rameters S m . 

In most experiments on single-molecule junctions the 
molecule is contacted by metal electrodes. Accordingly, 
we describe the left (L) / right (R) electrode by a contin- 
uum of non-interacting electronic states 



E 

keh/R 



Hur = } , £kC k c k 



(3) 



that are localized in the left/right lead. The coupling be- 
tween the molecule and the leads is given by 



#tun = J^ {V mk c\c m +h.c). 



(4) 



kEL,R;mEM 



where the coupling matrix elements V m k determine the so- 
called level- width functions 



rK,mn(e) = 27T ^ V mk V nkS(e - 6fc) 



(5) 



k£K 



with K G {L,r}. Throughout this work, we assume semi- 
infinite tight-binding chains as models for the leads with an 
internal hopping parameter 7 = 2 eV. The corresponding 
level-width functions are given by | 



r, 



K.mn 



(e) 



VK,mMK,n 

7 2 



V / 47 2 -(e- A i if ) 2 , 



(6) 
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where, similar to V m k, the parameters VK,m denote the 
coupling strength of state m to lead K. In addition, we 
assume a symmetric drop of the bias voltage <P at the con- 
tacts, i.e. the chemical potentials in the left and the right 
lead are given by /j.^ = e<P/2 and /ir = — e<P/2, respec- 
tively. Furthermore, we set the Fermi energy of the leads to 
£f = eV. The Hamiltonian of the overall system is thus 
given by the sum 



H — H e i 



H\ 



Hl + H R + H tl 



(7) 



For the NEGF approach that we introduce in Sec. 12.21 
it is expedient to remove the direct electronic-vibrational 
coupling terms in the Hamiltonian H by the small polaron 
transformation [3Q,l24l: 



H = e s He s = ^ e m c\ n c m + E fi««!,fl a 

m OL 

+ ^2(V mk X m c\c m + h.c.) + ^2 £kc\c k , 



(8) 



km 



with 



*J ^ / J \^ma.f "^ctj V^m,^Vn Vm)*ai 



ma 



X m = exp[iy2(\ ma /n a )P a ], 



(9) 



(10) 



and P a = —i(a a — a' a ). The effect of electronic- 
vibrational coupling is thus subsumed in a renormal- 
ization of the single-particle energies, e m = e m + 
(2S m — 1) J2 a (^rna/^a), the electron-electron interaction 
strengths, U mn = U mn - 2^ Q (A m aA„ Q //? Q ), and the 
molecule-lead coupling strengths that are dressed by the 
shift-operators X m . 

2.2 Nonequilibrium Green's Function Approach 
All single-particle observables, such as, e.g., the population 
of levels or the electrical current that is flowing through the 
junction, can be calculated from the single-particle Green's 
functions of the system. For the electronic degrees of free- 
dom of our present problem, they are given by 

G mn (r,T') = -i(T cCm (r)4(r')>H (11) 

= -<Cr c m (T)A m (r)cJ l (r , )Jft( T '))_ 

where the indices H/H indicate the Hamiltonian that is 
used to evaluate the respective expectation values. We em- 
ploy the ansatz l24ll25llT0ll2Tl: 



G mn (T,T') » G mn (T,r')(T c X m (T)Xl(T')) w , 



(12) 



to calculate these Green's functions with G mn (T,T') = 
— z(T c c m (T)c^(r'))-^ and T c the time-ordering operator 
on the Keldysh contour. The ansatz ( fT2l represents an ef- 
fective factorization of the single-particle Green's func- 
tions G mn into a product of the electronic Green's func- 
tions, G mn , and a correlation function of shift operators, 
(T c X m (T)X^(T')}ji-. This is justified, if the time scales 



for electronic processes and vibrational motion in the junc- 
tion are different, and is conceptually similar to the Born- 
Oppenheimer approximation OTI . 

We use an equation of motion technique to determine 
the single-particle Green's functions of the present prob- 
lem [24 25. 10.21 1. The equation of motion for the elec- 
tronic part of the Green's function reads 

(id T - e m )G m „(r, r') (-«?,-' - e„) = (13) 

6(T,T')5 mn (-id T > -e„) 

+£coul,mn(T) t') + 2_^ ^K,mn(T,T'), 
K£{h.R} 

where the self-energy contributions due to the coupling of 
the molecule to the left and the right leads are given by 

^ L /R, m „ (r,r') = (14) 

E K k Vn k gk(T, t0<Tc*»(t0*JLM>1? 
fceL/R 

with 5fc(r, t') the free Green's function associated with 
lead state k. The effect of electron-electron interactions is 
subsumed in the self-energy £coui,mn(T, t'). We describe 
the latter in terms of the elastic co-tunneling approxi- 
mation [32,33,10,2.1), where £cou\,mn(i~,T') is replaced 
by the self-energy £q>u1 mn( T > T ') describing electron- 
electron interactions in the isolated molecule (i.e. where 
V m k = 0), facilitating a nonperturbative description of 
electron-electron interactions. 

The correlation functions of the shift operators, 
(T c X m (T)X^(T'))jY, which represent the second part of 
the factorized Green's function G mn (cf. Eq. (fT2l). can be 
determined employing a second-order cumulant expansion 
in the dimensionless coupling parameters ^ff"- 11011211 

(T c X m (r)Xl(T / )) 7T = (15) 

v^ A ma A / . A ma A ,+A na A / 

„zZ QQ ' l -TTZ7T M r D cc,>( T > T !-■>■ ™n„Sl , m ^ D cc'(.T,T) 



with the momentum correlation functions 

Daa, = -i(T c P a (r)P a /(r'))ff. (16) 

These single-particle Green's functions can also be deter- 
mined evaluating their equations of motion 

{-d 2 r - n 2 a )D aa ,(T,T')(-d 2 T , - f2 2 a ,) = (17) 



AU £ Qj U L ' fV* 



1 



<5(t,t')<w(-^, - nl,)— + n el>aoi ,{T, T '). 

The corresponding self-energy matrix i7 e i jQ! Q' is evaluated 
up to second order in the molecule-lead coupling 01011211 



n, 



Amrv Ar 



d, Qtt '(T.O = -'L a a , 

J it r\j u _ r\>' 



(18) 



(S mn (T,T')G nm (T' ,t) + Z , „ m (r',r)G mn (T,r')). 

Since 77 e i, QQ ' depends on the electronic self-energies 
S mn = £ L , mn + ^R,mn and Green's functions G mn , 
the evaluation of Eqs. ( fT3l and ( fTTT i requires an iterative 
self-consistent solution scheme II10II211 . 
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2.3 Observables of interest To characterize elec- 
tron transport through a single-molecule junction, we an- 
alyze two different observables of interest as functions of 
the applied bias voltage <P: the average vibrational excita- 
tion and the electrical current flowing through the junction. 
The average level of excitation of vibrational mode a is 
given by 

(aia a ) H = (ala^ (19) 



E 



m 



■((c 



Sm)( 



•t, 



<5n)>H' 



where the latter terms represent the contribution from 
polaron-formation, which is associated with charging of 
the molecular bridge. It can be computed from the vibra- 
tional Green's function D aa according to ED 



1 



(a$ t a a )H*--Tm[D< a (t = 0)] 



(20) 



X 2 

V^ A ma 



in 



(11 

'E- 

m<n 



lm[G 



< 

mm 



-Im[G 



< 

mm 



--<■, 

(t = 0)] - 5, 



(t = 0)] - d m ) 

(lm[G<„(i = 0)]Im[G<„ 



(* = 0)] 



(lm[G< n (t = 0)]-S n ) 



function G mn and the 



where the electronic Green's 
Hartree-Fock-like factorization 

\ c m c mC n C n )jY ~ \ c m c m)\C n C n ) — \C, n C n )\C n C m ) (^-i) 

for m ^ n is used to represent the contribution from 
polaron-formation. 

The current through lead K, Ik, is determined by the 
number of electrons entering or leaving the lead in a given 
time interval (K <G {L,R}) 

Ik = -2ej f £ (c{c k ) w . (22) 

Here, the constant (— e) denotes the electron charge and 
the factor 2 accounts for spin-degeneracy. With the self- 
energies SuR.mn and Green's functions G mn (cf. Sec. 
12.21 ). the current can be calculated employing the Meir- 
Wingreen-like formula fJ34]|25][T0l 

Ik = (23) 



2e 



/ 



de 

2^ 



E 



^K.mn 



(e)G 



> ( 

■run \ 



y> ( f )C < 



(*)) 



Note that the NEGF approach is current conserving, i.e. 
Il = -h = I. 

2.4 Vibrations in thermal equilibrium In addition 
to the nonequilibrium treatment discussed above, we also 
employ a simpler 'equilibrium' description, where we ne- 
glect the vibrational self-energy, i.e. II e \, aa i(t,t') = 0. 
This is a commonly used approximation [30,24,25, 14 1 
and constitutes a powerful tool to elucidate and analyze 
nonequilibrium effects. Within this approximation, the vi- 
brational degrees of the junction are effectively confined to 



model 


ei 


H.,1 


fR,l 


7 


J7i 


An 


E1V1 


0.6 


0.1 


0.1 


3 


0.1 


0.06 


BAND 


0.15 


0.02 


0.02 


0.2 


0.1 


0.06 


REC 


0.6 


0.1 


0.03 


3 


0.1 


0.06 


RECBD 


0.15 


0.02 


0.006 


0.2 


0.1 


0.06 



Table 1 Model parameters for molecular junctions with 
a single electronic state (energy values are given in eV). 
The temperature in the leads is assumed to be T = 10 K 
throughout the article. Note that the parameters of all the 
model systems that we discuss in this review article fulfill 
the antiadiabatic condition Il/r.ii < Q\ and that they are 
in line with previous experimental H36I371 38 39 40l l2g]|4T1 
|5H and theoretical findings |42 !43 27 ,44 ,45 ,46) . 



the state they would acquire in thermal equilibrium at an 
effective temperature k^T. This means that, in contrast to 
a full nonequilibrium calculation, tunneling electrons may, 
indeed, excite or deexcite the vibrational degrees of free- 
dom of the junction but that such changes of the vibrational 
state decay instantly such that an electron that is tunnel- 
ing subsequently through the junction finds the vibrational 
mode in the same equilibrium state as before. 

3 Results We study the importance of resonant 
electron-hole pair creation processes for vibrationally cou- 
pled electron transport in single-molecule junctions in 
three steps of increasing complexity. First, in Sec. 13.11 we 
introduce the basic transport processes and mechanisms 
that occur in the resonant transport regime of a molecu- 
lar junction. To this end, we discuss the current-voltage 
characteristics of two generic model systems of molecular 
junctions that include a single and two electronic states. 
This provides the basis for our discussion in Sec. 13.21 
where we introduce electron-hole pair creation processes 
on the very same grounds as transport processes. Thereby, 
we analyze the vibrational excitation characteristics of 
the model systems that we already discussed in Sec. 13.11 
It is shown that electron-hole pair creation processes are 
crucial to understand the (current-induced) levels of vibra- 
tional excitation in a single-molecule junction, including 
the respective electrical transport properties. This is exem- 
plified in more detail in Sec. 13.31 where we discuss various 
transport phenomena that can be traced back to the influ- 
ence of electron-hole pair creation processes. Note that the 
validity of the results, which are presented in the following 
and based on the NEGF scheme outlined in Sec. 12.21 have 
been corroborated by comparison to Born-Markov master 
equation approaches 1101211[23l and a numerically exact 
wave-propagation scheme 11351 . 

3.1 Basic transport processes To discuss the ba- 
sic transport processes that occur in single-molecule junc- 
tions, we consider first a simple model for a molecular 
junction with a single electronic state that is coupled to 
a single vibrational mode (model El VI, see Tab. Q]for a 
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model 


ei 


£2 


Uia 


it,i 


H.,2 


^R,l 


^R,2 


7 


Qi 


Q 2 


An 


A21 


A12 


A22 


E2V1 


0.15 


0.8 





0.1 


0.1 


0.1 


0.1 


3 


0.1 


- 


0.06 


-0.06 


- 


- 


SPEC 


0.15 


0.8 





0.1 


0.1 


0.03 


0.03 


3 


0.1 


- 


0.06 


-0.06 


- 


- 


MSVE 


0.65 


-0.575 





0.1 


0.1 


0.03 


0.03 


2 


0.15 


0.2 


0.09 








0.12 


INT 


0.5123 


0.5126 


0.0248 


0.01 


0.01 


0.01 


-0.01 


0.25 


0.005 


- 


0.0248 


0.025 


- 


- 



Table 2 Model parameters for molecular junctions with two electronic states (energy values are given in eV). 
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Figure 1 Current- voltage characteristics of junction El VI 
(cf. Tab.Q]for the corresponding model parameters). 



detailed list of all model parameters). The corresponding 
current-voltage characteristic is shown in Fig. Q] where we 
distinguish three different scenarios: i) the electronic sce- 
nario (solid purple line), where we neglect the effect of 
electronic-vibrational coupling, ii) the vibronic scenario 
(solid black line), where we fully account for the coupling 
between the electronic state and the vibration, and Hi) the 
thermally equilibrated scenario (dashed black line), where 
electronic-vibrational is accounted for but the vibration is 
restricted to the state it acquires in thermal equilibrium at 
T = 10 K (which is effectively its ground state). 

Within the electronic scenario, the current-voltage 
characteristic displays a single step at e<P = 2e\, which 
is associated with the onset of resonant electronic trans- 
port processes where electrons are tunneling through the 
junction in two subsequent resonant tunneling events (see 
Fig. [2^). In contrast, the vibronic scenarios provide a num- 
ber of additional resonant transport channels that become 
active at different bias voltages 1141 . Some of them are 
readily active at the onset of the resonant transport regime 
at cP = 2ei and include both resonant excitation (Fig. [2)) 
and resonant deexcitation processes (Fig. [2};). Another set 
of resonant excitation processes (Fig. [2}i) becomes active 
at higher bias voltages, &P = 2(ei + nQ\) (n £ N), giv- 
ing rise to additional steps in the vibronic current- voltage 
characteristics. 

The height of the steps in the current-voltage character- 
istics is different for all three transport scenarios, reflecting 



the fact that the corresponding transport processes strongly 
depend on the electronic-vibrational coupling strength An 
and the vibrational excitation level of the junction. The 
respective probabilities can be expressed in terms of the 
so-called Franck-Condon matrix elements F mn , that is the 
transition matrix elements from the mth to the ?ith state 
of the vibrational mode upon change of the charge state 
of the molecule. If, for example, electronic transport pro- 
cesses (Fig. [2^) occur with a relative probability of 1 in 
the electronic scenario, they occur with a reduced prob- 
ability of _Foo ^ 0.7 in the thermally equilibrated trans- 
port scenario. This value becomes even smaller in the vi- 
bronic transport scenario, J^nPn^nn < Fqo, where p n 
represents the average population of the nth vibrational 
level. The first step in the current- voltage characteristics at 
e<P = 2ei is correlated with these probabilities, although a 
quantitative analysis is much more involved f 14|. The sub- 
sequent steps at e<P = 2(ei + nfi\) become successively 
smaller, which can be qualitatively understood by the re- 
duced Franck-Condon matrix elements of the respective 
transport processes that involve transitions from the vibra- 
tional ground state to the nth excited state of the vibration, 
F 0n = e~ A ii/ r2 i(Af"/^™)/n!. However, as the compar- 
ison of the dashed black and the solid black line shows, 
vibrational nonequilibrium effects significantly modify the 
step structure and lead, in general, to a suppression of the 
electrical current that is flowing through junction E1V1 

m. 

So far, we have considered transport processes for a 
model with a single electronic state. In many cases, how- 
ever, electron transport through single-molecule junctions 
is carried by a multitude of states. The resulting vibra- 
tional processes that may occur with respect to each of 
these states leads to a number of interesting (nonequilib- 
rium) phenomena (see, e.g., Refs. B10llllin~3lll4|[T7l ). As 
an example, we consider a junction with two electronic 
states (model E2V1, see Tab. [2] for a complete list of pa- 
rameters). The corresponding current- voltage characteris- 
tics are shown in Fig. [3] 

The current-voltage characteristic associated with the 
electronic scenario shows two steps at e<P = 2ei = 0.3 eV 
and e<P = 2e 2 = 1.6 eV, indicating the onset of electronic 
transport processes through states 1 and 2, respectively. 
Similarly, the thermally equilibrated transport shows three 
pronounced steps at e^ = 2|i = 0.228 eV, e<P = 2e 2 = 
1.528eV and e$ = 2(e 2 + U 12 ) = 1.672 eV, where the 
first of these steps is, just as in the electronic scenario, 
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(a) (b) (c) (d) 

LMR LMR LMR LMR 



r.r. 





Figure 2 Schematic representation of example processes for sequential tunneling in a molecular junction. Panel (a) shows 
sequential tunneling of an electron from the left lead to the right lead that involves two consecutive tunneling processes: 
from the left lead onto the molecular bridge and from the molecular bridge to the right lead. In panel (b)/(c) the latter 
of the two tunneling processes is accompanied by an excitation/deexcitation process, where due to electronic-vibrational 
coupling the vibrational mode is singly excited/deexcited. While the processes depicted by Panel (a) - (c) become active at 
the same bias voltage, that is for e<P re 2ei, resonant excitation processes like the one depicted by Panel (d) require higher 
bias voltages, e<P > 2ei + Q\. 




0. 0.5 1. 1.5 2. 2.5 

bias voltage O [V] 



Figure 3 Current- voltage characteristics of junction E2V1 
(cf. Tab. |2). 



associated with the onset of resonant electronic transport 
processes through the first electronic level. The other two 
steps originate from the onset of resonant electronic trans- 
port processes through the second electronic state. How- 
ever, due to electron-electron interactions, Un = 72 meV, 
transport through this state depends on wether the first 
level is occupied or unoccupied. This leads to a splitting 
of the corresponding step into two steps that are separated 
by eA<P = 2U\2- Apart from the splitting of steps due 
to electron-electron interactions, the current-voltage char- 
acteristics of the thermally equilibrated scenario can be 
understood on the same grounds as for a junction with a 
single electronic level. Vibrational side-steps are observed 
with respect to each of the three electronic resonances, 
that is at e@ = 2(ei + nQi), e<P = 2(e 2 + nfii) and 
e<£ = 2(e 2 + T7i2 + nQi) (n € N). They indicate, as be- 
fore, the onset of resonant excitation processes (see Figs. 
|2jl and 3k). This shows that the onset of vibrational pro- 




(b) 





Figure 4 Examples of sequential tunneling processes in- 
volving two electronic states. Panel (a) depicts a sequential 
tunneling process, which involves an excitation processes 
with respect to state 2. Thereby, it is assumed that state 1 is 
occupied, that is the tunneling electron requires an energy 
of 62 + U\2 + &i, which is increased by the charging en- 
ergy U12 due to interaction with the electron in state 1 . Pan- 
els (c) and (d) show corresponding deexcitation processes, 
where the lower-lying electronic state is unoccupied and 
occupied, respectively. 



cesses does not only depend on the position of the elec- 
tronic levels, given by e TO , but also on the electron-electron 
interaction strengths U mn . 

In contrast, the vibronic scenario shows a more com- 
plex behavior. For example, it exhibits only a single pro- 
nounced step at e<P = 2ei. The two steps associated with 
the onset of resonant transport through the second elec- 
tronic state are barely visible. This is due to the current- 
induced level of vibrational excitation, which is generated 
by inelastic processes with respect to the first state (Fig. 
|2}I). This leads, on one hand, to a reorganization of the 
steps heights and an overall reduction of the current level 
through state 1 (cf. the discussion of Fig.Q~]or see Ref. Ifl4l 
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Figure 5 Examples of electron-hole pair creation pro- 
cesses. Panel (a) shows a resonant electron-hole pair cre- 
ation process, where in two sequential tunneling events 
an electron tunnels from the left lead onto the molecular 
bridge and back again to the left lead. Thereby, the elec- 
tron takes up a quantum of vibrational energy (blue wiggly 
line). A pair creation process with respect to a higher-lying 
electronic state is shown in Panel (b). 



for a more detailed analysis of these phenomenona). On the 
other hand, it triggers resonant deexcitation processes with 
respect to the second electronic state (Figs.|5> and|4j;) even 
before it has entered the bias window H101I141 . As a result, 
the steps corresponding to transport through state 2 are ef- 
fectively broadened such that, instead of two pronounced 
steps, the onset of resonant transport through the second 
electronic state is associated with a number of small steps 
at e<P = 2(e 2 - nfii) and e<P = 2(e 2 + U 12 - nfii) ifTUl . 
This shows that vibrational nonequilibrium effects have a 
pronounced influence on the electrical transport propertie s 
of a molecular junction and that a more detailed and quan- 
titative understanding of these effects requires an analysis 
of the corresponding vibrational excitation characteristic 
(which will be given in Sec. 13.2b . 

3.2 Electron-hole pair creation processes In this 
section, we analyze the vibrational excitation characteris- 
tics of junctions El VI and E2V2. This requires to account 
for all processes, where the molecule exchanges energy 
with the electrodes, including inelastic transport processes 
(see Figs. [2] and |4|i as well as electron-hole pair creation 
processes (examples of which are depicted in Fig. |5). Pair 
creation and transport processes involve the same tunneling 
processes and, therefore, occur principally with the same 
probability. They are distinguished only by the fact that the 
tunneling electrons return to the same lead in the course of 
a pair creation event, while they transfer from one lead to 
the other in the course of a transport process. 

The excitation characteristics of junction E1V1 is de- 
picted by the solid black line in Fig. [6] It shows a quali- 
tatively different behavior with respect to the applied bias 
voltage <P than the corresponding current-voltage charac- 
teristics (cf. Sec. 13. Q , Although it exhibits steps at the 
same bias voltages, the height of these steps follows clearly 
another phenomenology. While the steps in the current- 
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Figure 6 Vibrational excitation characteristic of junction 
junction El VI (cf. Tab.Q}. 



voltage characteristic become successively smaller with an 
increasing bias voltage <P, they become successively larger 
in the vibrational excitation characteristic. Thus, in con- 
trast to the current-voltage characteristics, the level of vi- 
brational excitation does not saturate for high bias voltages 
l22l . Moreover, the step heights increase if the electronic - 
vibrational coupling strength is reduced. This can be in- 
ferred from the comparison of the vibrational excitation 
characteristic of junction El VI with the excitation char- 
acteristic of a very similar junction that differs from junc- 
tion El VI only by a reduced electronic- vibrational cou- 
pling strength An = 0.03 eV It is depicted by the solid red 
line in Fig. [6] This behavior is rather counter-intuitive, be- 
cause the transition matrix elements for the corresponding 
excitation processes (Fig. |2ji) become smaller for weaker 
electronic -vibrational coupling strengths An BOII47ll48l 

QUE). 

To understand the phenomenology of the vibrational 
excitation characteristics it is expedient to consider the na- 
ture of this observable. In contrast to the electrical current, 
which is determined by the number of transport processes 
in a given time interval, the level of vibrational excitation 
(in steady state) is determined by the ratio of the proba- 
bilities for excitation and deexcitation processes. At low 
bias voltages, e<i> < 2ei, vibrational excitation can only be 
generated in tunneling processes with respect to the right 
lead (cf. Fig. Wp)- Resonant deexcitation, however, is less 
restricted and can occur with respect to both leads {e.g. by 
the processes that are depicted in Figs. [2J; and [5^). This 
leads, roughly speaking, to a ratio between the probabilities 
for excitation and deexcitation of the vibrational mode of 
at least one to two. The corresponding level of vibrational 
excitation is, therefore, rather low and originates, to a large 
extent, from polaron formation (see Eqs. ([19} and (l20t). 
At higher bias voltages e^ > 2(ei + Q\), however, the 
imbalance between excitation and deexcitation processes 
becomes gradually smaller. This is, on one hand, due to an 
increasing number of excitation processes that become ac- 
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tive but, on the other hand, also due to a decreasing num- 
ber of deexcitation processes that were active at smaller 
bias voltages. For example, the electron-hole pair creation 
process depicted in Fig. [5^ is active at the onset of the 
resonant transport regime but is blocked for bias voltages 
e<£ > 2(ei + J?i). Thus, the probabilities for excitation and 
deexcitation of the vibrational mode become successively 
the same, resulting in an indefinite increase of vibrational 
excitation at higher bias voltages, where only transport pro- 
cesses occur l22l . Moreover, as the pair-creation processes 
that are blocked at e<£ > 2(ei + nfi\) are more important 
for weaker electronic-vibrational coupling, the correspond- 
ing increase/step in the vibrational excitation characteristic 
becomes larger for smaller values of An II14II221 . 

Similar as in Sec. 13. II we extend our considerations at 
this point to a junction with two electronic states, model 
E2V1. The vibrational excitation characteristic of this 
junction is shown in Fig. [7] As long as the second elec- 
tronic state is located far outside the bias window (for 
<P < 1 V), the level of vibrational excitation increases 
monotoneously due to resonant excitation processes with 
respect to the first electronic state (Fig.|2ji). At higher bias 
voltages, however, the vibrational energy thus generated 
facilitates resonant deexcitation processes with respect to 
the second electronic state, including transport (Figs. |4j> 
and Hj;) as well as electron-hole pair creation processes 
(Figs. \5fa). These processes significantly reduce the level 
of vibrational excitation before the second electronic state 
enters the bias window (i.e. leV< e<[> < 2?2)- The net 
effect is that the presence of the second electronic state 
results in a cooling of the junction at these bias voltages 

nun. 

This cooling extends to even higher bias voltages, if 
electron-electron interactions are present. This can be seen 
by comparison of the solid black line in Fig. [7] with the 
solid blue line, which shows the vibrational excitation 
characteristic of a junction that differs from junction E2V 1 
only by additional repulsive electron-electron interactions, 
U12 = 0.5 eV. Due to the more pronounced splitting of 
the resonances that are associated with the second elec- 
tronic state in this system, deexcitation processes due to 
both transport (Fig. |4j;) and electron-hole pair creation 
processes (Fig.|5j5) are shifted to higher energies and, con- 
sequently, become active at higher bias voltages. Thus, the 
additional electron-electron interactions lead to a slight 
increase of the level of vibrational excitation at low bias 
voltages (i.e. for <P < 2 V). At higher bias voltages, how- 
ever, they lead to a significant decrease of the vibrational 
excitation level (Coulomb Cooling) H141 . 

3.3 Transport phenomena due to electron-hole 
pair creation processes The results discussed in the 
preceding section, Sec. 13.21 show that a detailed under- 
standing of the vibrational excitation characteristic of a 
single-molecule contact can only be obtained if electron- 
hole pair creation processes are taken into account. Thus, 
pair creation processes also influence the corresponding 




1. 1.5 2. 
bias voltage <J> [V] 

Figure 7 Vibrational excitation characteristics for junc- 
tion E2V1 (solid black line, cf. Tab. 0. In addition, the 
vibrational excitation characteristics of a very similar sys- 
tem is shown, which, in contrast to junction E2V1, in- 
cludes additional repulsive electron-electron interactions, 
U 12 = 0.5 eV (solid blue line). 



electrical transport properties, as the efficiency of trans- 
port processes is strongly interrelated with the vibrational 
excitation levels (cf., e.g., the reduced current level in the 
resonant transport regime of junctions El VI and El V2 dis- 
cussed in Sec. 13. 11 1. In this section, we show that this influ- 
ence may be substantial and that it can lead to a number of 
interesting transport phenomena such as, for example, neg- 
ative differential resistance (Sec. 13.3. H i, rectification (Sec. 
13.3.2b . mode-selective vibrational excitation (Sec. l3.3.3l l or 
a pronounced temperature dependence of the current in the 
presence of destructive quantum interference effects (Sec. 

EH. 

3.3.1 Negative differential resistance In general, 
an electronic device exhibits negative differential resis- 
tance (NDR) when its differential conductance dI/d<P be- 
comes negative. This phenomenon is revelant for a num- 
ber of applications in electronics, for example, in analog- 
digital converters IJ491 or logic circuits I150JI511 . A prime 
example for a nanoelectronic device that exhibits NDR is 
the resonant tunneling diode, where the overlap of the (nar- 
row) conduction bands in the leads decreases with increas- 
ing bias voltage ^ [52 53]. The current-voltage character- 
istic of a corresponding model for a molecular junction 
(model BAND, cf. Tab. [T]) is depicted by the solid pur- 
ple line in Fig. [8] At low bias voltages, where the con- 
duction bands overlap with each other, the current level of 
the junction increases, in particular when resonant trans- 
port processes through the electronic state become active 
at e<S> ps 1t\. At higher bias voltage, however, the over- 
lap of the conduction bands continuously decreases such 
that electronic transport processes (Fig. |2^) can no longer 
connect an occupied state in one of the leads with an unoc- 
cupied state in the other lead. Thus, the current level of the 
junction decreases. 
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Figure 8 Current-voltage characteristics of a molecu- 
lar junction that exhibits negative differential resistance 
(model BAND, cf. Tab.[JJi. The inset shows the correspond- 
ing vibrational excitation characteristic. 



In the presence of electronic-vibrational coupling, the 
effect of the narrow conduction bands is significantly re- 
duced, as can be seen by the solid black line in Fig. [8] This 
is due to inelastic transport processes, which allow elec- 
trons to span a much wider range of energies in the course 
of a transport process, increasing effectively the overlap of 
the conduction bands. In addition, however, another much 
more narrow NDR feature occurs at e$ = 2(ei + fi\). 
It is in line with an increase of the vibrational excitation 
level of the junction (depicted in the inset of Fig. [8}. The 
thermally equilibrated scenario (dashed black line), how- 
ever, does not exhibit NDR at this bias voltage. On the con- 
trary, it displays an increase of the current level. Thus, we 
trace this NDR effect back to the suppression of electron- 
hole pair creation processes (Fig. |5^) and the associated 
increase in vibrational excitation rather than the opening 
of another inelastic channel (Fig.[2jl), although both occur 
at the same bias voltage. Note that this NDR effect is not 
related to the narrow conduction bands of this model sys- 
tem and that it may be more pronounced if the molecule 
is asymmetrically coupled to the leads (see, for example, 
Refs. l54ll55lfT4l). 

3.3.2 Vibrational rectification and spectroscopy 
in molecular junctions Diode-like behavior represents 
another important device characteristic ll37l . Such behav- 
ior is often associated with an asymmetric coupling of the 
molecule to the leads. But an asymmetry in the molecule- 
lead coupling is not sufficient l56l . as is demonstrated 
by the solid purple line in Fig. [9] It shows the electronic 
current-voltage characteristics of junction El VI with a re- 
duced coupling stength to the right lead (model REC, cf. 
Tab. [TJ. It is almost perfectly antisymmetric with respect 
to the polarity of the bias voltage, i.e. I(<P) = —/(—#). 
The origin of this behavior is twofold. First, transport pro- 
cesses involve the coupling to both leads and, therefore, 
cannot transfer an asymmetry in the coupling to the leads 



Figure 9 Current- voltage characteristics of the asymmet- 
ric junction REC (cf. Tab. [TJ- The inset shows the cor- 
responding vibrational excitation characteristic. The solid 
turquoise line refers to model RECBD (cf. Tab. [TJ, which 
is distinguished from junction REC in particular due to a 
much narrower band width 7. 



to the transport characteristics of a molecular contact. Sec- 
ond, the density of states for these processes does also not 
display a pronounced asymmetry with respect to the bias 
voltage <P. In particular in the limit ^r.i/Vl.i — > 0, where 
tunneling processes with respect to the right lead represent 
the bottleneck for transport, it is given by i~R,n(ei)/i / | i- 
If, however, in addition to an asymmetric molecule-lead 
coupling, the band width in the leads is reduced (from 7 = 
3eV to 7 = 0.2eV), the density of states -Tr,ii (ei)/^ 1 is 
much smaller in the resonant transport regime for negative 
bias voltages (e<P < 2?i) than it is for positive bias voltages 
(e<£ > 2?i). In that case, the corresponding current-voltage 
characteristic exhibits an asymmetry with respect to the po- 
larity of the applied bias voltage, which can be seen by 
the solid turquoise line in Fig. [9] which shows the current- 
voltage characteristic of such a junction (model RECBD, 

cf.Tab.EQlB 

Another mechanism that leads to a different current 
response at different polarities of the applied bias volt- 
age is electronic-vibrational coupling 01411211 . This can 
be seen by the solid black and dashed black lines in Fig. 
[9] which show the current-voltage characteristics of junc- 
tion REC considering the vibronic and the thermally equi- 
librated transport scenario, respectively. While the black 
line exhibits a pronounced asymmetry of the current level, 
the dashed black line shows an asymmetry in the current 
level only in the vicinity e^ ~ ±2ei. For the thermally 
equilibrated scenario, this behavior can be deduced again 
from the density of states for resonant transport processes 
in the right lead. For positive bias voltages, e<P > 2?i, it is 
given by Y,n F On r R,ii( e i)/^l.i ~ r R ,n(ei)/t/| a , where 



Note that for similar systems with ei ~ 0, the current re- 
sponse would be perfectly antisymmetric irrespective of an asym- 
metric molecule-lead coupling and the value of the band width 7. 
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Figure 10 Current- voltage characteristics of the asymmet- 
ric junction SPEC (cf. Tab. |2}. The inset shows the corre- 
sponding vibrational excitation characteristic. For compar- 
ison, the solid orange lines show results of the correspond- 
ing symmetric molecular junction E2V1 rescaled by a fac- 
tor of 1/5. 



the sum over the Franck-Condon matrix elements Fq u indi- 
cates that electronic (Fig.|2^) as well as resonant excitation 
processes (Fig. [J]}) become simultaneously active. For neg- 
ative bias voltages, however, only electronic processes are 
active at e<P — — 2ei, while excitation processes with re- 
spect to tunneling from the right lead become successively 
active at e<P = — 2(ei + nQ\). This leads to relative step 
heights in the current-voltage characteristics that are given 
by F 00 , Fqo + Foi, etc.. 

Although the same effects do also contribute to the 
asymmetric current response of the vibronic transport sce- 
nario, its behavior is strongly linked to the corresponding 
level of vibrational excitation (shown in the inset of Fig. [9]). 
As the excitation levels of the junction are much higher for 
negative bias voltages, the current is more strongly sup- 
pressed than for positive bias voltages. This relation be- 
tween high levels of vibrational excitation and a reduction 
of the current level has already been outlined in Sec. 13.11 
(and in Refs. 114112111 ). As a consequence, the current sup- 
pression or the rectification effect also extends in the vi- 
bronic scenario over a wider range of bias voltages than for 
the thermally equilibrated scenario. Thereby, the asymme- 
try in vibrational excitation is due to electron-hole pair cre- 
ation processes. In contrast to transport processes, they in- 
volve only one of the leads. Thus, in the limit ifc i/tt i — ► 
0, electron-hole pair creation processes with respect to the 
left lead (Fig. |5^) are much more important than pair cre- 
ation processes with respect to the right lead as well as 
excitation and deexcitation processes due to transport pro- 
cesses. Therefore, the vibrational mode is more efficiently 
cooled for positive than for negative bias voltages, since 
pair creation processes with respect to the left lead involve 
less vibrational quanta and, therefore, are more effective 
for positive than for negative bias voltages. 



Considering spectroscopic applications of single- 
molecule junctions, the above described vibrational rec- 
tification effect leads to important implications. For exam- 
ple, one may observe a different number of side-steps (or 
different relative step heights) in the current-voltage char- 
acteristic of an asymmetrically coupled molecular junction 
by changing the polarity of the applied bias voltage (cf. 
Fig. [9]l [5]. Furthermore, an asymmetry in the molecule- 
lead coupling may be useful to observe signals of states 
that are located further away from the Fermi level. Accord- 
ing to our discussion of Fig.[3l(Sec. l3. It . these signals may 
be blurred by the effect of resonant deexcitation processes 
ifTOl . In asymmetric junctions, however, where the effect 
of local cooling by electron-hole pair creation processes 
is different for different bias polarities, these processes 
can be less pronounced, at least for one polarity of the 
applied bias voltage 11231 . This is shown in Fig. [10] where 
the current- voltage characteristic of junction E2V1 is com- 
pared to a very similar junction that differs from junction 
E2V1 only by reduced coupling strengths to the right lead 
(model SPEC, cf. Tab. Q]). In contrast to the symmetric 
junction E2V1, junction SPEC exhibits, due to the effect 
of a more pronounced cooling by electron-hole pair cre- 
ation processes, two pronounced steps at e^ = 2ei and 
e<P = 2(?2 + U12) that correspond to the onset of reso- 
nant transport processes through the first and the second 
electronic state. 

3.3.3 Mode-selective vibrational excitation As 

our discussion of vibrational rectification effects in Sec. 
13.3.21 showed, the effect of cooling by electron-hole pair 
creation processes can be controlled by the external bias 
voltage <P if the molecule is asymmetrically coupled to 
the leads. Thereby, we focused on the electrical transport 
properties. In this section, we investigate junctions with 
multiple vibrational degrees of freedom and describe how 
the excitation levels of these modes can be selectively 
controlled exploiting the bias dependence of electron-hole 
pair creation processes B20ii2"T) . Such mode-selective vi- 
brational excitation may lead, for example, to applications 
in mode-selective chemistry ll57l [8l. where one seeks to 
break a specific (not necessarily the weakest) chemical 
bond in a molecule. 

To this end, we consider a generic model for a molecu- 
lar junction with two electronic states and two vibrational 
modes (model MSVE, see Tab. [2]). The two states are lo- 
cated above and below the Fermi level of the junction and 
may represent, for example, the HOMO and the LUMO 
level of the molecule. Moreover, the two states are coupled 
more strongly to the left than to the right lead. Such a cou- 
pling scenario is fairly common, in particular if one uses 
a scanning tunneling microscope to establish the contact 
to the molecule. The two modes are assumed to couple to 
one of the two states only, that is A mQ ~ S ma . This is 
an idealized situation, which allows to discuss the generic 
effect. Off-diagonal electronic-vibrational coupling dimin- 
ishes the effect, which, nevertheless, extends over a wide 
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Figure 11 Vibrational excitation characteristics of junc- 
tion MS VE (cf. Tab. 0. 



range of parameters (see Ref. lETl for a more detailed dis- 
cussion, including the effect of electron-electron interac- 
tions). 

The vibrational excitation characteristic of the two vi- 
brational modes is depicted in Fig.QT] In particular, mode 

1 shows significantly higher levels of excitation than mode 

2 for negative and vice versa for positive bias voltages, that 
is, the excitation levels of the two modes can be selectively 
controlled by the external bias voltage 02O1I211 . This can 
be explained by the same arguments that have been used to 
analyze the asymmetric excitation levels observed injunc- 
tion REC (cf. Sec. 13.3.2b . because the subsystem consist- 
ing of state 1 and mode 1 is effectively decoupled from 
the subsystem consisting of state 2 and mode 2. Thereby, 
it should be noted that the frequency of mode 2 is signifi- 
cantly larger than that of mode 1 . This means that for pos- 
itive bias voltages the (current-induced) vibrational energy 
is distributed in a way that is opposite to what statistical ar- 
guments would give in equilibrium systems. In more prac- 
tical terms, the stronger bond of the molecule, which is 
associated with mode 2, may eventually break in this junc- 
tion before the weaker bond. Note that such mode-selective 
vibrational excitation occurs also for other molecule-lead 
coupling scenarios I12TI . 

3.3.4 Temperature dependence of the current in 
the presence of destructive interference effects 
Electron-hole pair creation processes play also an impor- 
tant role if electron transport through a single-molecule 
junction is influenced by quantum interference effects |[T6l 
l4l [T71l . Quantum interference effects in molecular junc- 
tions are of interest not only from a fundamental point 
of view but also in the context of technological decive 
applications, including transistors [58], thermoelectric de- 
vices [59] or spin filters IJ601 . Interference effects occur 
in single-molecule junctions, e.g., when electron transport 
through the junctions is mediated by quasidegenerate elec- 
tronic states. In this section, we consider an example of 
such a molecular junction that includes a pair of quaside- 




#% 




*-* 



Figure 1 2 Panel (a) and (b): Graphical representation of a 
pair of quasidegenerate electronic states. Due to the differ- 
ent LoR symmetry of the two states, the tunneling ampli- 
tudes for (electronic) transport processes through these two 
states, which are depicted by purple wiggly lines, differ 
by a phase tt. Accordingly, they destructively interfer with 
each other, resulting in a suppression of the corresponding 
tunneling amplitude. 




20 40 60 80 100 

temperature in the leads [K] 

Figure 13 (Color online) Current-temperature characteris- 
tics of junction INT (cf. Tab. [2]). 



generate states: a bonding (symmetric) and an antibonding 
(antisymmetric) state (model INT, cf. Tab.|2]i. Note that the 
symmetry of the two states is reflected in the different sign 
of the molecule-lead coupling strength i/r.2- 

Due to the different symmetry of the two states, the 
tunneling amplitudes that are associated with electron tun- 
neling through the two states differ by a phase 7r and, there- 
fore, destructively interfer with each other (see Fig. \V2k 
andfTZb). This leads to a strong suppression of the tunnel 
current in this system ll6D . For example, at e<P — 0.2 eV, 
the electronic current level of junction INT is given by 
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~ 15pA0, while one obtains a current level of 0.4 fiA if 
interference effects are discarded in this system. 

As we have shown recently, electronic-vibrational cou- 
pling can strongly quench such interference effects |[T6l 
[171 . Thereby, it is decisive to note that the electronic states 
of a molecular junction couple differently to the vibrational 
degrees of freedom, even if they represent just symmetric 
and antisymmetric combinations of localized molecular or- 
bitals such as, for example, the states depicted in Fig. Q~2] 
Thus, due to the interaction with the vibrational degree of 
freedom, tunneling through one of the states becomes more 
favorable than through the other and, consequently, inter- 
ference effects are less pronounced. As a result, the current 
level of junction INT is significantly larger than without 
electronic-vibrational coupling. This can be seen in Fig. 
QjJ where the solid black and purple line depict the current- 
temperature characteristic of junction INT at e<P = 0.2 eV 
for the vibronic and electronic transport scenario, respec - 
tively. 

Thereby, the current level of the electronic scenario 
shows no temperature dependence because at bias voltages 
e<P 3> 2ei thermal broadening in the electrodes has no 
influence on the current level of the junction. In contrast, 
the vibronic current increases almost linearly with the tem- 
perature in the electrodes. This temperature dependence is 
due to the effect of electron-hole pair creation processes 
(Fig. [5]) [17]. As these processes involve only one of the 
electrodes, they are not suppressed by destructive interfer- 
ence effects. Therefore, the level of vibrational excitation 
is determined by electron-hole pair creation processes in 
this system rather than by inelastic transport processes that 
are suppressed by destructive interference effects. Similar 
as for a molecule that is adsorbed on a surface, pair cre- 
ation processes adapt the vibrational excitation level of the 
molecule to the temperature in the electrodes, leading to 
an almost linear increase of vibrational excitation, which 
is shown in the inset of Fig. [T3j as the temperature in the 
electrodes increases. As a result of the higher excitation 
levels, inelastic transport processes become more favor- 
able, leading to larger current levels in this system due 
to an enhanced quenching of destructive interference ef- 
fects. Thus, electron-hole pair creation processes facilitate 
a mechanism to control quantum interference effects in this 
junction by an external parameter, i.e. the temperature in 
the leads. This mechanism has recently been verified in a 
series of experiments on various single-molecule junctions 
by S. Ballmann et al. 0. 



4 Conclusions The theoretical studies 12511 10ll20lfT4l 
l22lfl6ll21ll23lfT7l of vibrationally coupled electron trans- 
port reviewed in this paper show that electron-hole pair 



2 Note that in this section, Sec. 13. 3.41 we refer to the electronic 
scenario as to the one without electronic-vibrational but including 
the polaron-shift of the electronic levels. Thus, the level of quasi- 
degeneracy between the electronic states in the vibronic and the 
electronic transport scenario is the same. 



creation processes play an important role in this nonequi- 
librium transport problem, although they are not directly 
contributing to the current that is flowing through the junc- 
tion. Similar to scenarios where a molecule is adsorbed on 
a surface, pair creation processes tend to adapt the vibra- 
tional excitation levels of the molecular bridge to the tem- 
perature in the electrodes. They are directly influencing the 
vibrational excitation levels of the junction and, thus, also 
the efficiency of transport processes. Although this repre- 
sents an indirect influence on the electrical transport prop- 
erties of a molecular junction it can, nevertheless, lead to 
a number of interesting transport phenomena such as, for 
example, negative differential resistance lfl4l and rectifica- 
tion 11411211 . In contrast to transport processes, electron - 
hole pair creation processes involve only one of the elec- 
trodes. As a consequence they transfer an asymmetry in 
the molecule-lead coupling to the corresponding transport 
characteristics. This is of relevance for spectroscopic appli- 
cations of single-molecule junctions l23l . mode-selective 
vibrational excitation [20,21] or in the presence of de- 
structive quantum interference effects [J4][T7] . Moreover, it 
was shown that the suppression of electron-hole pair cre- 
ation processes leads to higher levels of vibrational exci- 
tation in systems with weaker electronic-vibrational cou- 
pling lfl4l and that the absence of electron-hole pair cre- 
ation processes represents a link between the limit of high 
bias voltages and weak electronic-vibrational coupling (for 
<P > 2(ci + J7i)) 122]. 

Throughout the article, we have focused on resonant 
processes. It should be noted though, that, similar to res- 
onant and non-resonant transport processes [62. 25 1, off- 
resonant electron-hole pair creation processes are also of 
importance in molecular junctions. Aspects and implica- 
tions of these processes are discussed, for example, in Ref. 
ED . where it is shown that they play an important role in 
the non-resonant transport regime of a molecular junction 
but also at high bias voltages. 
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